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Abstract. We consider a class of geometrically frustrated Heisenberg spin systems 
which admit exact ground states. The systems consist of suitably coupled 
antiferromagnetic spin trimers with integer spin quantum numbers s and their ground 
state $ will be the product state of the local singlet ground states of the trimers. We 
provide linear equations for the inter-trimer coupling constants which are equivalent 
to $ being an eigcnstate of the corresponding Heisenberg Hamiltonian and sufficient 
conditions for $ being a ground state. The classical case s — > oo can be completely 
analyzed. For the quantum case we consider a couple of examples, where the critical 
values of the inter-trimer couplings are numerically determined. These examples 
include chains of corner sharing tetrahedra as well as certain spin tubes. $ is proven 
to be gapped in the case of trimer chains. This follows from a more general theorem 
on quantum chains with product ground states. 
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1. Introduction 

The effects of frustration in Heisenberg magnets liave attracted mucli interest over the 
past few decades [H [2] . Geometrically frustrated quantum antiferromagnets (AF) are an 
excellent play-ground for studying novel quantum many-body phenomena. We mention 
quantum spin-liquid phases, valence-bond crystal phases, order-by-disorder phenomena, 
lattice instabilities to name just a few. 

Moreover, in recent years there has been remarkable progress in synthesizing mag- 
netic materials [3]. Even exotic structures such as the star |4] or the maple-leaf lattices 
[5] have been synthesized. Hence, the investigation of exotic lattice structures being on 
the first glance purely academic might become relevant for experimental studies. 

The theoretical investigation of frustration effects in quantum spin anti- 
ferromagnets usually meets new difficulties; e. g. the quantum Monte Carlo method 
suffers from the sign problem for frustrated systems. Exact statements for interacting 
quantum-many body systems are rare and, therefore, new rigorous results are of con- 
siderable interest to improve their understanding. Moreover, solvable models may serve 
as test grounds for approximate methods. 

Starting with the seminal papers of Majumdar-Ghosh [6] and Shastry-Sutherland 
[7] , a famous class of spin systems with exact ground states has thoroughly been inves- 
tigated, see, for example, [8]. These systems consist of suitably coupled AF dimers 
such that the product state of the local 5 = dimer ground states remains the ground 
state of the total system. For this it is necessary that the coupling constants between 
the dimers satisfy certain linear equations [9] . A generalization to chains with trimerized 
ground states using SU{3) spins and bilinear-biquadratic Hamiltonians has recently been 
published [10]. Another approach using matrix product states as exact ground states 
[TT] has also been mainly applied to anisotropic spin systems. 

In this paper we will rather consider isotropic (Heisenberg) spin systems and ex- 
tend the basic idea underlying the Majumdar-Ghosh-Shastry-Sutherland ground states 
to systems with AF trimers as building blocks. These systems will be geometrically frus- 
trated. If the individual spin quantum number s is integer, the ground state of a uniform 
(or almost uniform) AF trimer will be again a non-degenerate singlet state, i. e. having 
S = 0. The product state $ of these local ground states will be an eigenstate of the 
Heisenberg Hamiltonian of N trimers if and only if the coupling constants satisfy cer- 
tain linear constraints completely analogous to the dimer case. Moreover, if the coupling 
between the dimers is not too strong, $ will be the ground state of the total system, 
called the trimerized ground state (TGS). Of course, the question arises what precisely 
is meant by "not too strong"? In general, the precise domain of systems admitting TGS 
can only be numerically investigated for given examples and lattice structures. In a 



Coupled trimers 



3 



sense, we thus invert the usual strategy to numerically find ground states for given spin 
systems. We define a certain state $ and numerically calculate the coupling constants 
for systems which have $ as their ground state. However, we have also derived some 
general rigorous statements on systems with trimerized ground states. These statements 
mainly concern trimerized ground states in the classical limit s — )■ cxo, which can be com- 
pletely understood, and some sufficient conditions for trimerized ground states in the 
quantum case. Moreover, we prove that in the case of trimer chains with trimerized 
ground state, $ will be gapped for an extended domain of coupling constants. We stress 
that our investigation is restricted to the case of integer s, since for s being half integer 
the ground state of the AF trimer will be degenerate, see, e. g. [12]. We mention that 
trimerized states of the kind explained above have also been considered as approximate 
ground states of certain modified 2D kagome lattices [T3] . 

The paper is organized as follows. In section |2] we summarize the basic definitions 
and main results. Section [3] is devoted to a couple of examples, starting from a single 
trimer, followed by a pair of trimers and then passing to 1-dimensional chains of trimers. 
The latter examples cover systems like chains of corner-sharing tetrahedra [H] and spin 
tubes [15], which currently attract a lot of attention in the literature, albeit, as a rule, 
not under the aspect of exact ground states. The study of the examples led us to 
some conjectures about the domain of the coupling constants for systems admitting 
trimerized ground states and its dependence on s and A^. In order to clearly distinguish 
these conjectures from the mentioned rigorous statements we have presented the latter 
in a separate section HI We decided to give detailed proofs of these statements only in 
those cases where they markedly differ from the analogous proofs for the dimer case in 



2. Basic definitions and summary of main results 

We consider systems of 3A^ spins with one and the same individual integer spin quantum 
number s = 1,2,3, .. . which are grouped into N fixed triples ("trimers"). To indicate 
this grouping the spins will be denoted by indices /i = (z, 6) where z = 1, . . . , is the 
trimer index and 6 = 0, 1, 2 distinguishes between the three spins belonging to the same 
trimer. Further we consider Heisenberg Hamiltonians 



where = (s}* , sfj, , sfj, ) denotes the /i-th spin observable and J the 3A^ x 3A^- matrix 
of real exchange parameters or coupling constants J^,^ satifying 




(1) 




(2) 



i,j <5,<: 



J. 



IMIJ. — 



0, J, 



Jy^ for all /i, = 1, . . . , 3A^ . 



(3) 
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All operators act on a (2s + 1) ^^-dimensional Hilbert space n = (^IZ^ n^. If the spin 
quantum number s is fixed, we may identify a spin system with its matrix J. Note that 
due to the matrix notation each scalar product ■ occurs twice in the Hamiltonian 
([1]). In order to comply with the usual notation we have therefore introduced the factor 
i in some examples of section |3l 

For any trimer with index i let [iO, il, i2] denote the ground state of the AF trimer 
-^0 = Ylse ' '^«<:; Aj > whlch Is uulque up to a phase factor. If the trimer index i 
is irrelevant, it will be simply denoted by [0, 1, 2]. For general s this state can be written 
in terms of the Wigner-3j-symbol as 

[ s s s \ 

[0,1,2]= |mo,mi, -mo - mi) , (4) 

^ — ' \ mn mi — mn — mi / 

mo,m\=—s \ / 

using the eigenbasis |m), m = — s, . . . , s of s/t and the corresponding product bases. 
The state [i0,zl,i2] will remain the ground state of the trimer i even when its Hamil- 
tonian is suitably disturbed, see section 13.11 Let Cf denote the set of all ( J^o, Jn, Ji2) 
where this is the case. 

The ground state of a system of unconnected AF trimers satisfying ( J^o, Jn, Ji2) G 
Cf for z = 1, . . . , A^ is the product state 

TV 

$^ = (g)[^0,2l,^2] , (5) 

called the trimerized state; it has the total spin quantum number S = 0. A system J is 
said to admit trimerized ground states (TGS), or to have the TGS property, iff is a 
ground state of -^(J), i- e. iff 

($^1^(1)$^) < (6) 

for all G H with ||\E'|| = 1. Let C|, denote the set of all spin systems J with the TGS 
property. If the quantum number s is understood, we suppress it and write simply $ 
and 

For J G C$ it is necessary that $ will be an eigenstate of H{S). This turns out to 
be true if and only if the inter-trimer coupling constants fulfil the relations 

JiOJO + Jiljl = JiOJl + Jil,jO (7) 

■JiOJO + <^tlj2 = Ji0,j2 + JiljO (8) 

■JiOJO + "^12^1 = -JiOJl + "^12,^0 (9) 

JiOJO + ^i2j2 = JiOJ2 + Ji2J0 (10) 

for all 1 < i < j < A^. The corresponding eigenvalue 

E = -S{S + 1) J^O,il + J^0,^2 + ^^1,^2^ • (11) 
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is independent of the inter-trimer coupling. Since (17|)- (fT0|) is a system of four independent 
linear equations, the set of all real, symmetric 3A^ x 3 A^- matrices satisfying fl7|)- f|T0|) and 
Jfiii = for all = 1, . . . , 3A^ will be a linear space of dimension 3A^ + 5(^) = y(5A^+l), 
denoted by J7$. The set C$ of TGS systems will form a convex cone embedded in the 
linear space ^7$, since the condition ([6]) is invariant under positive linear combinations 
of J's , see also [9]. At the boundary of C$, the trimerized ground state $ will become 
degenerate, i. e. there will exist "competing" ground states which have a lower energy 

o 

than (fTTIl if J crosses the boundary of Cj,. Sometimes we will also use the symbol C$ 
in order to denote the open convex cone of spin systems J where $ is a non-degenerate 
ground state. 

The conditions (j7l)-( fT0l) still include interesting spin structures such as corner- 
sharing tetrahedra and spin tubes to be considered in section 13.31 In some of these 
chains additional symmetries arise which allow the description by an equivalent ladder 
model of composite spins, see 13.3.11 and 13.3.21 However, the TGS property is indepen- 
dent of this additional symmetry as shown by the example of a certain spin tube in 
section 13.3.21 

It is clear that the choice of dimensionless numbers for the J^^ in the examples 
implies the introduction of appropriate units for energy and temperature in order to 
apply the Heisenberg model to real systems. In this sense, temperature becomes a di- 
mensionless quantity in the thermodynamic calculations of section 13.21 which illustrate 
some physical consequences of the existence of TGS ground states. In most cases the 
domain of TGS systems, i. e. the shape of the cone C|, can only be determined numer- 
ically. Exceptions are the single trimer case where C|, can be calculated analytically, 
see section [3711 and the two-trimer case with s = 1, see section [3l2l In all examples 
which are considered in section [3] there is some evidence that these cones shrink with 
increasing s. Note that in the classical limit will be degenerate, since the necessary 
conditions for the J^^ are stronger in this case, see section I4l2l In contrast, for the chains 
of trimers considered in section 13. 3[ the dependence of C|, on is only weak. More- 
over, we found that C|, is slightly expanding if grows. This is an indication that the 
limit of C|, for iV — )■ oo is "non-degenerate", in the sense that it does not converge to a 
lower- dimensional domain and that the trimerized ground state is gapped. We will rigor- 
ously prove these properties in section l44l Our class of models thus supports Haldane's 
conjecture [16] that integer spin chains possess a unique ground state and an energy 
gap between the ground state and the excited states. The famous AKLT model [17] is 
also an s = 1 spin chain with a unique gapped ground state but its anti-ferromagnetic 
Heisenberg Hamiltonian is modified by a biquadratic term. For the s = 1/2 Heisenberg 
spin chain the absence of an energy gap has been proven in the 'Lieb-Schultz-Mattis 
theorem' [T8] . 
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In the classical case s — )> oo an analogous definition of TGS systems is possible, 
see section HI In this case $ consists of all spin configurations with a mutual angle 
of 120° between spin vectors of the same trimer. For classical TGS systems the inter- 
trimer coupling must be necessarily uniform and hence can be described by a symmetric 
X N— matrix G. We have the result that a classical system has the TGS property if 
and only if G is positive semi- definite. Thus the classical case is completely understood. 

3. Examples 

For readers less interested in the mathematical details of our rigorous analysis that 
will be presented in more detail below in section H] we first present some examples. In 
particular, the general statements listed above can be used to discuss certain specific 
chain-like models, such as chains of corner sharing tetrahedra as well as various spin 
tubes. We begin with some more elementary examples. 

3.1. One Trimer 

We consider three spins with integer spin quantum number s and Heisenberg 
Hamiltonian 

3 

H = ^ J/^uSf, ■ = J1S2 ■ S3 + J2S3 ■ si + J3S1 ■ S2 , (12) 
/i,v=i 

where we have relabeled the coupling constants in order to keep the following formulas 
readable. Let [1,2,3] denote the unique state with vanishing total spin, S = 0. It is an 
eigenstate of (1121) . since the eigenspaces of S"^ are invariant under H. For certain values 
of Ji, J2, J3, [1,2,3] is even the ground state of H, e. g. for Ji = J2 = J3 = 1- These 
values of Ji, J2, J3 form an closed convex cone in the 3-dimensional (Ji, J2, J3)-space. 

o 

At the boundary of C the ground state of H becomes degenerate. Recall that denotes 
the open subset of where [1,2,3] is the non-degenerate ground state. The form of 
the cone can be calculated using computer-algebraic software and the well-established 
assumption that the competing state has the quantum number S = 1. It is given by the 
following inequalities 

(Ji,J2,J3) eC^ ^ 

1 9 

-S{S + l)(Ji + J2 + .h){JlJ2 + J2-h + ^3^1) < (1 + 7^S{S + 1)) Jl J2J3 

and Ji, J2, J3 > . (13) 

The intersection of C'^, s = 1, . . . , 10 with the plane Ji -|- J2 + J3 = 3 is depicted in figure 
[TJ By Taylor expansion of (fT3|) one can show that, for increasing values of s, the cones 
approach a circular form centered at the half line Ji = J2 = J3 > with a cone angle 

of 




Figure 1. Representation of the neighborhood of Ji = J2 = = 1 in which the state 
[1, 2, 3] remains the ground state of ([T^ where s = 1, . . . , 10 starting with the outermost 
curve. The values of Ji, J2,^3 are restricted to the plane Ji + J2 + J3 = 3. The 
curves have been calculated according to (jl3p and have been confirmed numerically for 
s ~ 1, . . . , 5. The point at the center of the figure with coordinates (1,1,1) corresponds 
to the classical limit. 

Hence the cones will shrink and approach their classical limit Ji = J2 = J3 > for 
s — )■ 00, see section m 

3.2. Two Trimers 

Next we consider two trimers, i. e. six spins grouped into two triples with indices (1, 2, 3) 
and (4, 5, 6) and Heisenberg Hamiltonian 

6 

fi,i'=i 

We ask whether $ = [1, 2, 3] ® [4, 5, 6] will be a ground state of H. $ is then called the 
trimerized ground state. First, we note that $ need not be an eigenstate of H unless 
the Jfj,i, do not satisfy equations fl7|)-f lT0|) . These equations will be derived in section HI 
They can be expressed in the following way: Let {X, Y) be a pair of spins belonging to 
the first trimer and (x, y) another pair belonging to the second one. Then the two sums 
Jxx + Jvy and Jxy + Jyx niust be equal. This is a kind of balance condition completely 
analogous to the corresponding condition in the case of two dimers, see [9]. It is satisfied 
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for all combinations of spin pairs if and only if $ is an eigenstate of H. Actually, only 
four conditions of the above form have to be postulated, since the other will follow then, 
see section HI The coupling constants within the same trimer are not constrained, hence 
we are left with a 3 + 3 + 5 = 11-dimensional linear space of independent coupling 
constants J^j,, which will be called J^. It is independent of s. 

Again, the linear space will contain a closed convex cone of values J^^, such 
that $ will be a ground state of the corresponding Hamiltonian f lT5|) . The example 
of two unconnected trimers with the respective non-degenerate ground states [1,2,3] 
and [4, 5, 6] shows that C** is not empty. But it seems hopeless to analytically calculate 
except for s = 1. Thus we will here present some mixture of numerical, computer- 
algebraic and semi-analytical results. Neglecting a positive overall factor in f|T5|) we still 
have a 10- dimensional manifold of possible J^j,- values which is difficult to visualize. We 
will hence confine ourselves to some two-dimensional subspace of defined by (see also 
figure [2]) 

Ji2 = Ji3 = -hs = J45 = J46 = J56 = ^ and (16) 

_ T — T— ^T— ^T— T— ^^'^ 

Jl4 = — -'26 = = 2' "^15 — 2' "^25 — -'SS — — ^ — 5 (1' J 

where obviously the second equation fulfils the general conditions (17|)- (fT0|) . Note that 
the coupling strength between spins indicated in figure |2] equals twice the values of the 
Jfj_i, due to our definition of the Hamiltonian (fT5|) . 



The set T'^ of points with coordinates (a, b) such that $ will be a ground state of the 
corresponding Hamiltonian is again a convex set. It is represented in figure [5] for the 
values s = 1,2,3,4,5. The case s = 1 is just within the practical limits of computer- 
algebraic methods. The boundary of consists of pieces of 6 intersecting curves given 
by equations of the form qij{a, b) = where the qij are polynomials in the variables a, b 
with integer coefficients. The three simplest cases are 

goo = 4 - 4a - + - 46 + 5ab - a% - 6^ - ab'^ + b^ , (18) 
go2 = 8 + 28a + 28a^ + 8a^ -8b- Uab - 5a^b - 2b^ - bab^ + 2b^ , (19) 

and 

gio = 128 - 128a - TAa^ + 85a^ - 2a^ - lla^ + 2a^ - 1286 + 226a6 - 29a=^6 

- 68a^6 + 25a^6 - 2a^6 - 746^ - 29a6^ + 66a^b^ - 8a%^ - 2a%^ + 85b^ 

- 68ab^ - 8a%^ + Aa%^ - 2b^ + 25a6^ - 2a%'^ - llb^ - 2ab^ + 2b^. 

(20) 

The other polynomials are too complicated to be reproduced here. Note that the Hamil- 
tonian (fT5|) with the coupling constants (fT6 |[T7|) commutes with 5^^ and 6^23- Correspond- 
ingly, the indices i,j of the polynomials qij refer to the quantum numbers S = i and 
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Figure 2. Two coupled trimers satisfying (fT6)) . (|17p . The intra-trimer coupling (black 
lines) is set to one. The inter trimer coupling strengths are a (red line), —a (green 
line), b (blue line) and b — a (magenta line). 

5*23 = j of the corresponding competing states, see figure [3l 

In order to illustrate the physical implications of the presence of TGS ground states 
for a relatively simple example we have calculated the (dimensionless) zero-field mag- 
netic susceptibility x = (1^)1^-0 ~ ^ function of temperature T for coupling 
constants along the line b = OA a and s = 1, see figure HI This line crosses the boundary 
of the TGS domain at the value ai = 0.869506904299778. For a < ai the ground 
state $ has = 0; for a > ai the competing ground states have 5 = 1. Hence there 
is a transition from x{T) vanishing exponentially at T = to divergence of the form 
x{T) ~ 1^- The factor | is simply the mean value of for the three ground states 
with 5" = 1 and M = —1,0, 1. For a = ai the factor is |, corresponding to the mean 
value of for the four ground states with 5" = 1 and M = —1, 0, 1 and S = M = 0. 
This transition is qualitatively the same for gapped infinite TGS chains, see section l44l 
although the factor | would have to be replaced by the mean value of of a continuum 
of competing ground states. 

In the two-trimer example the sets are shrinking when s increases. We generally 
conjecture that D C^' for s < s' but could not prove this rigorously. For 6 = and 
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b 

1.5 r 




Figure 3. Representation of the set T*, s = 1, of points with coordinates {a,b) 
such that the singlet product state $ will be a ground state of the Hamiltonian (fTS)) 
corresponding to (|16ll7p . The boundary of consists of pieces of 6 intersecting 
smooth curves qij which have been calculated by computer-algebraic means. Note 
that qii contributes twice to the boundary. The indices i,j of qij refer to the quantum 
numbers S = i and 5*23 = j of the corresponding competing states. The dotted line 
defined by 6 = 0.4 a will be used for thermodynamic calculations, see figure HI 



a < we find numerically that the critical value acrit which lies at the boundary of 
has the form 

acrit = ^ for 5 = 1,2,3,4,5. (21) 

s + 1 

This can be confirmed semi-analytically to hold for all s by calculating the competing 
state which becomes an additional ground state if a assumes the value flcrit- We defer 
this calculation to the Appendix. 




Figure 4. The product of temperature T and zero-field susceptibility x as a function of 
T for the coupling constants in (|17p varying from a = 0.8 (lowermost curve) to a = 0.9 
(uppermost curve) and h = 0.4 a. The units are chosen such that Tx = (M^) becomes 
dimensionless. At the value ^ 0.869506904299778 (black curve) the line & = 0.4 a 
crosses the boundary of the TGS domain and consequently there is a transition from 
to the finite value 1/2 resp. 2/3 of lim7-_^o Tx due to the finite magnetization of the 
competing state with S = 1, see figure [H 



3.3. Chains 

Next we consider examples of chains formed of trimers which are coupled in a balanced 
way, i. e. satisfying (!7|)- (fT0|) . such that the singlet product state 

N 

<i> = (^[i0,il,i2] (22) 

1=1 

becomes an eigenstate of the corresponding Hamiltonian. The coupling within the 
trimers is always chosen as 

JiO,il = Jil,iO = JiQ,i2 = Ji2,i0 = Jil,i2 = Ji2,il = " • (^3) 

The coupling constants between the trimers Js^e = JiS, {i+iy are chosen to depend linearly 
on one or two parameters and we will investigate the convex domain of these parameters 
for which $ will be a ground state of the corresponding Hamiltonian. Throughout this 
section we adopt periodic boundary conditions, i. e. + 1 = 1. 



3.3.1. Chains of corner sharing tetrahedra 



The choice of the inter-trimer interaction matrix 

^ 6 0^ 
(J^J = ^ I 6 

b J 



(24) 
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Figure 5. Representation of tlie set of points with coordinates (a, b) such that the 
singlet product state $ will be a ground state of the Hamiltonian (fT5|) corresponding 
to ()16|17p . The different curves defining the boundary of T" refer to the values of 
s = 1,2,3,4,5 starting from the outermost curve. The curve for s = 1 has been 
calculated analytically, see figure [3l the other curves for s = 2,3,4,5 have been 
determined numerically. 



leads to a chain of comer sharing tetrahedra, see figure [6l Similar systems have been 
widely considered theoretically as well as experimentally, see e. g. [H]. However, most 
of these studies are focussed on spin-half systems. As it is typical for chains of corner 
sharing tetrahedra the chains considered in this section have an additional symmetry, 
namely that certain composite spin squares (in figure [6] these are 6^235 §56' ^^c- ) com- 
mute with the Hamiltonian. Hence we have a simpler sawtooth chain composed of spins 
and composite spins, see figure [6l While in general the composite spins may have spin 
quantum number Scomp = 0, . . . , 2s, in the singlet product ground state $ the composite 
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(2,3) (5,6) (8,9) (11,12) 



6 



12 



Figure 6. Example of a chain consisting of corner sharing tetrahedra (left) and an 
equivalent effective chain of composite spins (right). 



spins have the same spin quantum number s as the individual spins. 

We have numerically determined the critical values feniin ^-^d femax where $ ceases 
to be the unique ground state with energy Eq = —^s{s + 1). These critical values 
depend on the number of spins 3A^ and the spin quantum number s. The results are 
contained in table [H While the dependence on the length of the chain is weak, again 
an increase in the spin quantum number leads to a significant smaller parameter region 
where $ is the ground state. From the numerical data we can detect the competing 
states which become ground states for b < fe^iin ^ > Viax- The competing state for 
negative, i.e. ferromagnetic, 6 is a ferrimagnetic state (the total spin of the chain 5* is 
finite but less than 3iVs) and the composite spins have the spin quantum number s — 1. 
Hence for s = 1 the spin quantum number of the composite spins in the competing state 
is zero. Then the spins along the base line of the effective sawtooth chain build a simple 
ferromagnetic s = 1 chain which is decoupled from the composite spins. The energy of 
the competing state is Eq = N{—2 + b) and its total spin is 5" = A^. As a result there is 
a prominent transition for s = 1 at b^i^ = — 1. For s > 1 and feromagnetic b there is 
no simple competing state, since the effective chain is a mixed-spin sawtooth chain [19] 
(e.g. a mixed spin-one spin- two sawtooth chain in case of s = 2). 

The competing state for positive, i.e. antiferromagnetic, 6 is a state with the total 
spin of the chain S = and the composite spins have the spin quantum number s — 1. 
Again for s = 1 the spin quantum number of the composite spins in the competing 
state is zero, and, as a result, the spins along the base line of the sawtooth build an 
antiferromagnetic s = 1 Haldane chain which is decoupled from the composite spins. 
The energy of the competing state is Eq = N{—2 + beo{N)). Setting Eq = Eq yields 
^max = ~ eo{N) ■ N ^ oo the energy of the antiferromagnetic chain is given by 
eo(oo,l) = -1.40148403897, see [20]. Hence we get ioi N ^ oo the critical value 
^max = — l/eo(oo,l) = 0.713529353. This result together with b^^^ = —1 indicates 
that the domain of TGS systems slightly expands with growing A^. Hence we expect a 
non-degenerate TGS domain even for N oo where $ will be a gapped ground state as 
we will prove in section 14.41 For s > 1 the effective chain is again a mixed-spin sawtooth 
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Table 1. Critical values of the coupling constant b corresponding to the chains of 
figure ini 



3N 


s 


^min 


^niax 


12 


1 


-1.0 


0.667 


12 


2 


-0.57 


0.434 


18 


1 


-1.0 


0.696 


24 


1 


-1.0 


0.706 


30 


1 


-1.0 


0.710 



chain and one can find values for femax for sliort cliains only. 

Finally, we have numerically calculated the energy gap A^v of the chain of N corner 
sharing tetrahedra as a function of the coupling constant b for = 2, 3, 4, 5, see figure 
[71 These results confirm the rigorous bound An > A2 for all = 3, 4, . . ., which will 
be derived in section 14.41 corollary [2] and shows that the TGS ground state is gapped 
for trimer chains of this kind. 




{Js,e) = o I —a —a 
—a —a 



3.3.2. Spin tubes Several specific choices of the inter-trimer interaction matrix between 
the trimers correspond to so-called three-leg or triangular spin tubes. Such spin tubes 
have been widely considered in the literature, see e. g. [T5] . 
Spin tube I: The choice of the inter-trimer interaction matrix 

\ 

(25) 

/ 

leads to a spin tube as shown figure [HI For the special kind of systems ([25]) it turns 
out that the chains have an additional symmetry, namely that certain composite spin 
squares (in figure [8] these are 6^23, S\^, etc. ) commute with the Hamiltonian. Hence one 
can consider a simpler ladder model composed of spins and composite spins, see figure 
[H This ladder is frustrated, since the upper and the lower leg exchange bonds have 
different sign. In the singlet product state $ the composite spins on the lower leg have 
the same spin quantum number s as the individual spins. Hence we have an effective 
spin-s ladder, where the antiferromagnetically coupled rungs are in a local singlet state. 
Note that a similar exact singlet product state was also found for a spin-half ladder with 
one ferromagnetic leg and one antiferromagnetic leg |21j . 



We have numerically determined the critical values amin and Omax where $ ceases to 
be the unique ground state with energy Eq = —^s{s + 1). These critical values depend 
on the number of spins 3N and the spin quantum number s. The results are contained 
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b 



-1.0 -0.5 0.0 0.5 1.0 

Figure 7. The numerically calculated energy gap Ais[{b) of the chain of N corner 
sharing tetrahedra as a function of the coupling constant b for TV = 2, 3, 4, 5 and 5 = 1. 
The lower- most blue curve corresponds to N = 2; the next curves correspond to = 3 
(red), = 4 (black) and N = 5 (green). According to corollary 2 of section we 
have the rigorous bound An > A2 for all A'^ = 3,4, . . ., which is confirmed by these 
examples. The curves Ajv(6) show only small variations for A'^ = 3, 4, 5. They seem to 
grow monotonically with A^ except for b > 0.64 where the A^ = 3 and the A^ = 4 curves 
intersect. The parabolic form of the gap functions in the neighborhood of 6 = can be 
understood by virtue of 1st and 2nd order perturbation theory and considering local 
S = 1 excitations for 6 = satisfying Ajv(O) ~ 1. The 1st order corrections vanish 
due to the balanced form of the inter-trimer Hamiltonian; the 2nd order contributions 
have a negative sign since the local excitations can be treated like ground states in the 
Hilbert space Note that the perturbed local excitations need not be the competing 
ground states; this explains the kinks in the gap functions. 




Figure 8. Example of a spin tube I (left) 
composite spins (right). 



and 



an equivalent effective ladder of 
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Table 2. Critical values of the coupling constant a corresponding to the chains of 
figure IHl 



3A^ 


s 


'^min 




12 


1 


-0.319 


0.418 


18 


1 


-0.330 


0.419 


12 


2 


-0.210 


0.270 



in table [2j In accordance with our general conjectures we find that the critical interval 
[amin, Omax] shriuks when passing from s = 1 to s = 2 and slightly expands from = 12 
to iV = 18. 



Based on the numerical data we have analyzed the competing states which become 
ground states for a < Omin and a > Omax- In the competing state for both cases the 
composite spins have the spin quantum number s + 1. As a result, the lower leg of the 
effective model carrying larger spins determines the magnetic ordering of the systems. 
While for positive a within both legs the spin-spin correlations are ferromagnetic in the 
competing state, one has a competing state with antiferromagnetic spin-spin correlations 
within both legs for negative a. Due to the antiferromagnetic rung coupling the total 
spin of the system is S" = for a < a^nin, whereas the competing state is ferrimagnetic 
with S* = for a > flmax- 



Spin tubes II: The choice of the inter-trimer interaction matrix 
/ d 0\ 



d 
d 



(26) 



leads to another spin tube, see figure [91 For this special kind of coupling it turns out 
that the chains have no additional symmetry (i.e. no composite spins are conserved), if 
both parameters b and d are non-zero. Nevertheless, $ will be the unique ground state 
for a convex neighborhood of the point b = d = 0. Our numerical results are contained 
in figure [TOl 

First, we notice that in the limits b = or d = some of the bonds are missing, 
and the model can be transformed to the chain of corner-sharing tetrahedra. As a 
trivial result the competing states and the corresponding transition points are the same 
as discussed in section 13.3.11 Consequently, the competing state at (and also in the 
vicinity of) b = 0, d = —1 and also d = 0, b = —1 is ferrimagnetic. In all the other areas 
we find numerically that the competing state is a non-trivial singlet state. However, its 
spin-spin correlations depend strongly on the position on the transition line shown in 
figure [TOl 




0.8 
0.6 
0.4 
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-0.2 
-0.4 
-0.6 
-0.8 
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-3N=12: s=l - 
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3N=12: s=2 - 
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3N=18: s=l 


o 

1 1 


1 1 1 1 



-1 -0.8 -0.6 -0.4 -0.2 



0.2 
d 



0.4 0.6 0.8 1 



Figure 10. Representation of the set of points with coordinates {b,d) such that $ 
win be a ground state of Hamiltonian corresponding to spin tube II shown in figure |9l 



As for the previous examples, it is obvious from figure [TO] that there is only a very 
weak dependence on the size of the system. Moreover, we observe an inclusion C T^, 
which is compatible with the conjecture that the domain in the coordinate space (6, d) 
where $ will be a ground state shrinks with increasing spin quantum number s. 



4. Rigorous results 
^.i. Definitions 

We recall the general definitions given in section |2j Analogous definitions hold for the 
classical case: Here the spin observables s^' are unit vectors, if(JJ)'^' is the Hamiltonian 
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function, defined on the 3A^-fold Cartesian product of unit spheres 

3N 



V 



X si) , 



(27) 



and C P is the set of all spin configurations satisfying 

Sio + Sii + Si2 = o for alH = 1, . . . , . (28) 

Note that ^'^^ as well as are invariant under rotations. J is said to have the classical 
TGS property iff the minimum of H{S) is assumed for all s G ^'^K In this case we write 



4-2. Necessary conditions for TGS systems 

Whereas a complete characterization of C|, seems to be possible only for small N and 
s or for the classical case s = oo, one can prove a number of partial results, either 
necessary or sufficient conditions for J G 

We have already mentioned the following result which gives a necessary condition 
for J G Q: 

Theorem 1 ^ is an eigenstate of H{J) iff 



JiOJO + Jiljl 


JiOJl + JilJO 


(29) 


■'' + Jil,j2 = 


Ji0,j2 + Jil,jO 


(30) 


JiO,jO + "^12^1 = 


JiO,jl + "^12^0 


(31) 


JiO,jO + ^t2j2 = 


Ji0,j2 + Ji2,j0 


(32) 



for alii < j = 2, . . . , N . Moreover, let "DSq denote the space of all real 3 x 3-matrices with 
vanishing row and column sums, and Jij the 3 x 3-matrix with entries Jiejs, e, 5 = 0, 1, 2. 
Then the above four equations - are equivalent to the statement that Jij is 
orthogonal to the space VSq w. r. t. the inner product {A,B) = Tr (A^B). 

Proof: The second part of the theorem follows, since the equations (!29|) - (!32l) say that 
Jij is orthogonal to the four matrices 



/ 1 

-1 




-1 o\ 

1 




V 



/ 1 



-1 



-1 o\ 



1 



v 



1 

-1 



/ 1 
-1 
^00 

which span VSq- 

To prove the first part of the theorem we rewrite the Hamiltonian (!27l) in the form 

H{^) = 5m ■ 5m (34) 



(33) 



i^j £,5=0 
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N 

(35) 

j=i 

= E¥.' (36) 



where the distribution of the terms of the second sum in ( l35l) to the terms H is 
arbitrary. We have H = Hij ® 1 such that acts on =1-Li® Tij and 1 on 
the remaining factors. Recall that the trimerized state has the form 

N 

<l> = (g)[20,^l,z2] , (37) 

i=l 

where [z0,zl,z2] denotes the AF trimer ground state in "Hj = Tiio (g) Tin (g) 'Hi2- The 
following lemma can be proven completely analogous to lemma 2 in [9] 

Lemma 1 ^ is an eigenstate of H{J) iff' [i0,il,i2] ® [j0,jl,j2] is an eigenstate of Hij 
for alii < j = 2,...N. 

In view of this lemma we only need to consider the case of = 2 trimers with 
indices i < j in the remaining part of the proof. We set = [iO, il, i2] and rewrite the 
indices according to 

(,0) = 1, (tl) = 2, (^2) = 3, (jO) = 4, (jl) = 5, (j2) = 6 . (38) 

Since all summands in 

= (015^0) = T.(<P\is? + 5? + ^3^)"<f>) (39) 

i=l 

are non-negative, we conclude ||(si*'* + Sg"* + Sa"*) 0|P = 0, i. e. {si^ + + Sj^^) = 
for i = 1,2, 3. Further, (si + S2 + S3) ■ s^(0 ® 0) = for i/ = 4, 5, 6. Hence, for arbitrary 
^ G Hij, the matrix D with entries 

= (^'|s^-s,+3(0®0)) /i,i/ = l,2,3, (40) 
has vanishing row and column sums, i. e. D E "DSq. ® is an eigenstate of H[J) 
iff ^(0 ® 0) = with H = Y.l.M=i J^^u Sfj_ ■ 8^+3- This in turn is equivalent to 
(^|4(0 ® 0)) = for all ^ e or {J,D) = Y.l,u=i Jf^'^D^, = 0, i. e. J is or- 
thogonal to all matrices in VSq which can be written in the form (l40l) . 
It remains to show that there exist enough \1/ G Hij such that the matrices of the 
form (HOj) constitute a basis of T>Sq. Note that for all s = 1, 2, 3, . . . has a non- 
vanishing scalar product with the basis vector e = |1, —1,0). Choose = e e and 
consider the corresponding matrix (HOl) D^^^ with entries oj^J = ■ Si/+3(0 <S) 0)) = 

(e|s^0) (e|s^0). First, we conclude that {e\sj!'^(f)) = {e\sj^^(j)) = since Sj}'' and s^f-* 
change the spin quantum number S'^gs which is for e and 0. Second, Sg e = 0, 



hence D^^^ = di 



( 1 -1 o\ 

-1 1 

y y 



di 7^ 0. The remaining matrices of the basis in (l33l) 
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are similarly obtained by choosing \1/ = |1, —1, 0, 1, 0, —1), = |1, 0, —1, 1, —1, 0), and 
= |1, 0, —1, 1, 0, — l). This concludes the proof of theorem [H ■ 

Since (!29|) - (!32l) is a system of four linearly independent equations, the set of all real, 
symmetric 3N x 3iV-matrices satisfying ( |29l) - (j32l) and J^^ = for all yU = 1, . . . , 3A^ will 
be a linear space of dimension 3N + 5(^) = ^(SiV+l), denoted by J7$. The set C$ of 
TGS systems will form a closed convex cone embedded in the linear space J7$, see [9]. 
If $ is an eigenstate of if(J) it is straightforward to calculate the corresponding 
eigenvalues, since ■ = for i 7^ j: 

Corollary 1 If ^ is an eigenstate of -ff (J) then 

^(J)$ = -s{s + 1) J,o,il + J^0,i2 + Jil,^^ $ ■ (41) 

In the classical case we have similar but stronger results: The conditions fl29|) -fl32l) 
can be strengthened to a uniform coupling condition: 

Theorem 2 // J G then the coupling constants do not depend on S,e, i. e. 

JiS,ie = Ji>0 (42) 

and 

Ji&,je = ^ij (43) 

for all 5, e = 0, 1, 2 and i < j = 2, . . . , N . 

Consequently, we will denote by = J7|' the linear space of all real, symmetric, 
3N X 3A^-matrices J with vanishing diagonals and satisfying (l42l) (except Ji > 0) and 

Proof of theorem O The ground states of classical Heisenberg systems satisfy 

^ J^^s^ = K^s^, = 1,...,3A^ , (44) 

V 

see eq. (16) in [22]. This equation results from the condition that the energy H = 
' assumes a minimum, subject to the constraints ■ = 1. Here 
the K^, fi = 1,...,3N appear as the Lagrange parameters corresponding to these 
constraints. We choose fi = {i,0) and rewrite ( l44l) in the form 

(45) 

It is clear, by definition of classical TGS systems, that the contributions from different 
trimers with index j ^ i in (145|) can be rotated independently. These rotated contri- 
butions cannot be compensated by variations of unless Jioje^je vanishes for all 
j 7^ i. Choosing s^g^ = —sfi = \ and sjg^ = yields Jioji = Jio,j2- Similar arguments 
apply to the other equations which say that the coupling between different trimers must 
be uniform. 

To prove uniform coupling within the trimers we reconsider ( l45l) in the form Jio^nSn + 
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Jio,i2Si2 = HioSiQ. The special choice s-|^ = —s^^ = | and = again yields 
Jio,ii = Jio,i2 and analogously for the other equations. 

The previous considerations show that for classical trimerized ground states the Hamil- 
tonian assumes the value E = —3 Xlili "^j- If '^^^ of Jj would be negative, say Ji < 0, 
one could lower the energy by choosing siq = Su = Si2. Hence all Jj > and the proof 
is complete. ■ 



4-3. Systems close to unconnected trimers 

o 

Let J denote the matrix of an unconnected TGS system, i. e. 

o o o o 

( Jio.ii, Jio,j25 Jii,j2) for alH = 1, . . . , where denote the local trimerized ground 

o o o 

states. All other matrix elements Ji^je with i ^ j vanish. Of course, JIgC|, and hence 
the next lowest energy eigenvalue Ei satisfies 

El = Eo + e = +e, e > . (46) 

o 

By continuity arguments, a small neighborhood of J still consists of TGS systems. 
We want to derive a more quantitative result and consider an inter-trimer 3A^ x 3A^ 
symmetric coupling matrix A 7^ which has to be "small" in a certain sense. As a 
measure of "smallness" of A we will use |5min| where (5min denotes the lowest eigenvalue 
of the matrix A. Note that Tr A = 0, hence (5min < and the highest eigenvalue (5max 
of A satisfies < 5max ^ (3A^ — l)|5min|- It is clear that the size of the neighborhood of 

o o 

J depends on the energy gap of iJ(JI) which explains the e in the numerator of ( l47ll : 

o o 

Proposition 1 Let JI he an unconnected TGS system and JJ =J +A, A G such that 

l"^™"! ^ 3Ns{s + l) ' ^^'^^ 
where 5min denotes the lowest eigenvalue of A. Then J G CJ. 

Although the proof of proposition [1] is largely analogous to that of proposition 3 in [9], 
we will give it here for sake of convenience. The s-dependence of the bound in fjTTj) 
supports the conjecture that the cones C| shrink with increasing s. 
Proof of proposition [U Let \E' be any normalized state satisfying \E' ± $. It follows that 



Further, 



(^lif(J)l^) > El . (48) 
(^|if(A)|$)= ^A^,, ($|s^-s,|^) = TtAS (49) 



fJ,,V 



3N 

\ ^ /iTr I „2 I 



> d^in TtS = 5^in 2^(^|s^|^) = 3Ar<5^ins(s + 1) (50) 
>-e, (51) 
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where the last inequahty follows from (H7|) and (Jmin < 0. It follows that 

{^H{^)\^) = (*|i/(J)|*) + (*|iJ(A)|*) (52) 
>E,-e = Eo, (53) 
hence $ will be a ground state of iJ(J). This concludes the proof of proposition [H I 

An important special case of proposition [1] is the case of an unconnected 

o o o 

homogeneous TGS system, i. e. Jio,ii=Jio,i2=Jii,i2= Aj > for all i = 1, . . . , TV. In 
this case 

e = 2\ = 2mm{Xi\i = l,...,N} . (54) 

One of the simplest potential TGS systems J(e), see figure [TTl shows an interesting 
effect: For given s and sufficiently small e it is a TGS system by virtue of proposition [1] 
But if e > is fixed and s increases, it eventually looses the TGS property. Otherwise 
we would get a contradiction since JI(e) ^ by theorem [2] and the (normalized) ground 
state energy must converge for s — ?■ oo towards its classical value as a consequence of 
the Berezin/Lieb inequality [23] 

{s + ifE^J^,^ < E^^ < s'Et^ . (55) 



4-.4- The TGS chain 

A TGS chain consists of copies of trimers such that the intra-trimer and inter-trimer 
coupling is invariant under one- dimensional translations and $ is a ground state with 
energy Eq{N). A first question is whether for a given coupling the system remains a 
TGS chain for all A^ G N, A^ > 2. If this is the case, one may ask if the difference 
An between the next-lowest eigenvalue Ei{N) and Eo{N) has a positive lower bound 
independent of A^. In this case one says that $ is a gapped ground state. Here we ignore 
further questions concerning the limit of A^r for N ^ oo and confine ourselves to the 
existence of a gap for TGS chains. 

We will prove our result in a slightly more general context. Correspondingly, some 
of the general definitions of the paper are abolished in the following theorem. We 
consider a Hamiltonian H = Xlili where the Hi live in Hilbert spaces "Hj ® T-Li+i 
and the total Hamiltonian H in Ti = <S)iLi Hilbert spaces Hi are copies of 

one finite-dimensional Hilbert space. These and the following definitions are to be 
understood in the sense of cyclic boundary conditions A^ -|- 1 = A^. Moreover, if T is 
the unitary translation operator in shifting the tensor factors cyclically and hence 
satisfying = 1, we will assume T HiT* = Hi^i and hence [T, H] = 0. 

Theorem 3 Let $j G Tii be normalized and ^i (g) be the unique ground state of 
Hi with eigenvalue Eq and the next-lowest eigenvalue being Ei = Eq + 6, 6 > 0. Then 
= <S>i=i be the unique ground state of H with eigenvalue Eq = N Eq and the 

next-lowest eigenvalue satisfies Ei > Eq + 25 . 
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Figure 11. This system cannot be a TGS system for fixed e > and arbitrary s, since 
its classical limit is not TGS. The classical ground state corresponding to an energy 
i?o = — 3 — %- is indicated by small arrows. 



Proof: The first claim follows immediately by 

N N 

H<l>= Hi^ = ^ $1 ® . . . ® Hi{^i ® ®---^N 

i=l i=l 

N 

= Eo$ = A^^o'^' = ^o'^' (56) 

i=l 

and NEq being an obvious lower bound of H. 

Let G H be the eigenvector of H belonging to the next-lowest eigenvalues Ei > Eq. 
We may assume that 

* ± $ and = e^"^, a = 27ik/N, k eZ . (57) 

Our aim is to show Ei > Eq + 26. Let /i = 0, 1, 2, . . . denote the eigenbasis of 
Hi in "Hj ® 'Hi+i such that |0) = $i ® \fJ',K) denotes a corresponding product 

basis in T-L, where K stands for some multi-index of quantum numbers. Moreover, we 
consider the reduced density operator in "Hj ® Hi+i defined by the partial trace 

(/x|iy4|j.)=5^(/i,/i|M/)(M/|z/,A'). (58) 

K 
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Then we conclude 

N 



El = = (59) 

1=1 

N 

= 5^Tr(i7,H^4) (60) 

i=l 

= J2T^T{E,\ix){f^\W^) (61) 

i,ti 

= J^EMWM (62) 
= J2 i Eo{Om\0) + Yl E,{fi\W^\fi)] (63) 



^l=l,2,... 

>Y,(Eo{0\Wm + {Eo + 5) J2 if^l^M] ■ (64) 

i \ 11=1,2,... J 

Lemma 2 (0|VI^4|0) < 1 - ^ • 

For the proof of the lemma we use an arbitrary orthonormal basis |?t,), ?t, = 0, 1, 2, . . . 
in Hi such that |0) = and T operates as a cyclic shift operator in the corresponding 
product basis in "H. Hence (0|iy|,|0) will be rewritten as (0, 0|W^|,|0, 0). In this notation 
we have $ = |0, 0, . . . , 0). Due to translational symmetry the term (165|) does not depend 
on i, hence we may take i = 1 in what follows. We conclude 

{0,Om\0,0) = Y,mO,0,K)\' (66) 

K 

and 

1 = W4= Yl I(*|0,0,n3,n4,...)|'+ Yl mni,n2,...)\' 

n3,n4,... ni,n2v 
= So + Si . (67) 

The first sum sq in ( !67j) runs through all sequences 0,0,77,3,77,4, . . . excluding the value 
773 = 774 = . . . = 0, since = 0. Equivalently, we will say that it runs through 

all states ip = |0, 0, 773, 774, . . .) G Bq. The second sum si in ( l67l) runs through all 
sequences 771,712, .. . except those with 77i = 712 = 0, or, equivalently, through all states 
= |77i,7Z2, . . .) G Bi. Thus the total sum in (|67l) runs through an orthonormal basis 
B = BoUBiof W = {tfj e n\{i/j\^) = 0}. 

We consider on B the equivalence relation ~ '?/^2 ^1 = "^2; a G Z, and denote by 
A = B/^ the corresponding set of equivalence classes or "orbits". Due to (!57l) all states 
ip in the same orbit A yield the same value 

tA = |(*|V^)|' = |(*|r'^^)|2, aGZ. (68) 

For each orbit A G A let Nx = |A| denote its length. For most orbits we have Nx = N, 
but in general Nx will be a divisor of N. For example, if iV = 6 and |1, 2, 3, 1, 2, 3) G A 
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then A^A = 3. We define N\ = \\r\Bk\., /c = 0,l, and obtain the following equations: 

iVA = iVf ) + iv{') , (69) 

^0 = Y.t,Nf\ (70) 

AeA 

AeA 

If N\ = N any basis vector ip = |?7,i, 77-2, tt-s, 77.4, .. .) G A has exactly mutually 
orthogonal translations. Note that at least one Uj, 1 < j < N must be non-zero since 
7^ $ = |0, 0, . . . , 0). Hence at least two translations of ip belong to Bi, namely those 
where j is shifted to 1 or 2. It follows that A^{^^ > 2 and hence N^^ < N—2. Similarly, in 
the general case of 1 < Nx<N we also have N^^^ > 2 and hence N^^^ < Nx-2 < N-2. 
In the case A^a = 1, that is, ip = \n,n, . . . ,n), n > we have A^^"*^^ = 1 and A^^"'' = 0. 
This case has to be treated separately. We write A G Ai iff A^a = 1 and A G A> iff 
A^A > 1 and conclude 

^0= J^tAiVf <(A^-2) J]tA, (72) 

A6A> AgA> 

Sl= 5Z ^A + J] tAiVf > 2 ^ tA , (T3) 
AeAi AgA> AeA> 

which for ^a6A> ^a > implies 

- > ^ • (74) 
So ~ - 2 ^ ' 

If ^A6A> ^A = then So = and (l65l) follows immediately. From ( FM|) we infer 
and 

s^>^, (77) 
which concludes the proof of the lemma. I 

To complete the proof of theorem [3] we consider 

l = TTW^ = {0\Wi\0)+ J2 {l^\WM = so + si (78) 

tM=l,2,... 

and rewrite flM|) as 

N N 

Ei>Yl (^0 ^0 + (^0 + S) si) = Y^iEo + S si) (79) 

i=l 1=1 

= NEo + N6si>NEo + N6^ = Eo + 26, (80) 
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where we have used ( I77I) which is equivalent to (1651) . ■ 
The generahzation of theorem [3] (ci = 1) to square {d = 2) and cubic {d = 3) lattices is 
obvious but will not be considered here. We only note that in this case the energy gap 
is bounded from below by 2"^ 5 for = 1, 2, 3. 

In order to apply theorem [3] to trimer chains we will take 

1=1 i=l S,e=0 ^ 

+ JseSi^S ■ , (81) 

where the Jse satisfy the conditions of theorem [1] and the Hilbert spaces Hi, i = 1, . . . , N 
are chosen appropriately. Of course, $i = [io,ii,i2]- 

The 3 X 3-matrix J contains five independent real numbers and may thus be considered 

o 

as a vector of M^. We will fix the values of the intra-trimer coupling Jse such that the 
open convex set 

T = {J\Hi has the unique ground state [iQ,ii,i2]} C (82) 

is non-empty. Hence the energy gap oi Hi, Ei — Eq = S{J) varies over T but remains 
positive there. Then theorem [3] shows that H^^^ remains a TGS chain for all values 
J e T and all N e N. Moreover, the energy gap An{J) = Ei{N) - Eo{N) satisfies 
AAr(J) > 2 6{J), hence $ is a gapped ground state in this case. Note that, due to the 
cyclic boundary conditions, we have H^"^^ = 2Hi, hence A2(J) = 26{J). We summarize: 

Corollary 2 Let H^'^'> be the Hamiltonian of a trimer chain according to [8l\) and the 
intra-trimer coupling J he chosen such that [8^] is non-empty. Then H'^^^ will he a 
TGS chain for all J E T and all N = 2,3, . . . and its unique ground state $ possesses 
an energy gap satisfying An{J) > A2(J). 

4-5. The classical case 

In the classical case it is possible to completely characterize all TGS systems. Recall 
that eij = eji denotes the uniform interaction strength between two trimers and Jj that 
within the trimers according to theorem [2l For any J G J7|' we define an x A^- matrix 
G(J) with entries 

Gii = Ji for alH = 1, . . . , AT , (83) 
Gij = eij for alH ^ j = 1, . . . , iV . (84) 
Then we have the following result: 

Theorem 4 Let J G JT"!', then J G C|f iff is positive semi-definite. 
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Recall that G(J) > iff the principal minors det(Gjj)jj=i^...^„, > for n = 1, . . . ,N. 
Hence for classical spin systems the TGS property can be checked by testing N inequal- 
ities. 

This result is also relevant for quantum spin systems, since we have the following: 
Proposition 2 C C|, for all s = 1,2,3, .. .. 

Again, the proofs of theorem |4] and proposition [2] are analogous to those given in sections 
5.4 and 5.7 of [9]. 

Appendix: Another eigenstate for the trimer pair 

We reconsider the system of two trimers in section 13.21 with the special coupling 



and the remaining coupling constants vanishing. We want to calculate a competing 
eigenstate which gives a lower energy than the trimerized state $ for a < acrit- 

As usual, we denote the composite spin of a subsystem by subscripts, e. g. S^fgs = 
(si + S2 + Ss)^ with eigenvalues S'i23(-S'i23 + 1). Recall that it is possible to construct 
orthonormal bases Bt for the Hilbert space belonging to a spin system by means of 
"coupling schemes" T. For example, the coupling scheme 1 —t- 12 —t- 123 yields the 

common eigenbase of the composite spin squares SI2, §123 and the 3-component of the 

(s) 

total spin 5*^23 • happen that some state belongs to different Sy's. In this 

case will be a common eigenvector of all composite spin squares corresponding to the 
different coupling schemes T. For example, the 5 = ground state [1, 2, 3] of a uniform 
AF trimer belongs to both coupling schemes 1 — j- 12 — 123 and 1 — )■ 13 — j- 123 and 
hence has the good quantum numbers Su = S13 = s. (Otherwise it cannot couple with, 
say, 5*3 = s to give S'123 = 0). Note that a vector \E' will belong to different Bts if it is 
already uniquely determined by a proper subset of the quantum numbers corresponding 
to some coupling scheme Ti. If T2 is any other scheme containing the same quantum 
numbers we have necessarily \1/ G Bt2- 

We will apply these considerations to the coupling scheme Ti = (2 — )• 23 — )■ 123, 5 — )■ 
56 — )■ 456 —J- 123456) and the quantum numbers S'23 = S^q = s + 1, S'123 = 1, S'123456 = 
S = 0. Obviously, it follows that 5*456 must have the value 1 and the vector \& is uniquely 
determined by these quantum numbers. The same subsystems 23, 56, 123, 123456 also 
occur in the coupling scheme T2 = (5 ^ 56,2 ^ 23 ^ 123 1234 123456). Hence 
\1/ G Bti n Bt2 and thus \1/ must have the quantum number S'1234 = s + 1. Analogously, 
>S'i456 = s + 1. It follows that is an eigenstate of the Heisenberg Hamiltonian 



J12 — Ji3 — J23 — J45 — Jm — <^56 — 2' 

a 

JlA = — J26 = —^36 = —^25 = —^35 = 77 



(85) 



a < , 



(86) 
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_ ( q2 _ C2 _ c2^ _ f C2 _ C2 _ c2\ 
I £1234 £123 £4 ) I £l456 £456 £1 ) 

= 2(-si • S4 + S2 ■ S5 + S2 ■ S6 + S3 ■ S5 + S3 ■ se) (87) 

with eigenvalue Ei = — 4(s + 1). If H{a) denotes the Hamiltonian according to f fT6|) . 
( IT7|) and 6 = 0, we conclude that 

H{a) = (gLs - Ms + 1) + §456 - 3s(s + 1)) - a //i , (88) 

and hence {^\H{a)\'^) = 4a(s + 1) +4 -6s(s + 1). This equals {^\H{a)\^) = -6s(s + l) 
for a = Ocrit = "7^4 which confirms fl2T]) . Note that the argument is not completely 
rigorous, since we could not exclude other competing ground states than for arbitrary 
s. However, it is in agreement with our numerical data presented above. 
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